Phase behaviour of additive binary mixtures in the limit of infinite asymmetry 
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We provide an exact mapping between the density functional of a binary mixture and that of the 
effective one-component fluid in the limit of infinite asymmetry. The fluid of parallel hard cubes is 
thus mapped onto that of parallel adhesive hard cubes. Its phase behaviour reveals that demixing of a 
very asymmetric mixture can only occur between a solvent-rich fluid and a permeated large particle 
solid or between two large particle solids with different packing fractions. Comparing with hard 
spheres mixtures we conclude that the phase behaviour of very asymmetric hard-particle mixtures 
can be determined from that of the large component interacting via an adhesive-like potential. 
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Since Biben and Hansen showed [Gj that a binary mix- 
ture of hard spheres (HS) in Rogers- Young approxima- 
tion has a spinodal instability for diameter ratios smaller 
than 0.2, it is generally accepted that a demixing tran- 
sition occurs in this fluid if the asymmetry between the 
two component sizes is large enough. This fact has been 
confirmed with other theories |2| , simulations || , and ex- 
periments in colloidal suspensions Q . The question is no 
more whether such a mixture demixes but which is its 
phase behaviour. Some approximate theories |^] suggest 
that the phase rich in large spheres is a crystal instead 
of a fluid, in agreement with what is observed in the ex- 
periments [Q . The depletion potential between the large 
spheres due to the small ones has been determined both 
perturbatively |6| and in simulations jTj]. If the inter- 
action of the small spheres is replaced by this potential 
the resulting effective fluid can be simulated (7[j|]. In 
the unstable region of the binary mixture large clusters 
of depleted spheres are formed, which evolve very slowly 
0. The phase behavior of this effective fluid for diam- 
eter ratios smaller than 0.1 shows a large coexis- 
tence region between a diluted fluid and a highly packed 
solid, and even an isostructural solid-solid transition, if 
the small spheres packing fraction and the diameter ratio 
are low enough. 

The picture emerging from these studies is that the 
phase behaviour of a very asymmetric binary mixture 
can be understood in terms of the phase behaviour of a 
fluid of HS interacting via a depletion potential. This po- 
tential has a very narrow (essentially the diameter of the 
small spheres) and deep well [0J|], and this type of po- 
tential have been proven to have no liquid phase JlO|,[ll| : 
there is only fluid-solid coexistence, or — if it is narrower 
and deeper — in addition there may appear an isostruc- 
tural solid-solid transition. 

In this letter we will explore this connection between 
the phase behaviour of an additive binary mixture and 
that of the effective depleted one-component fluid by 
explicitly showing that the Hclmholtz free-energy func- 
tional of the mixture reduces to that of the correspond- 
ing one-component adhesive fluid in the limit of infinite 
asymmetry. That in this limit the solvent induces an ad- 



hesive interaction between the large particles is already a 
known result 12 RJ; but we are going a step further by 



mapping the functional of the mixture onto that of the 
corresponding adhesive fluid, and hence their respective 
phase behaviours. 

The procedure is as follows. As usual in these systems 
we fix the chemical potential of the small particles, ps 
(semigrand ensemble), and therefore the system is de- 
scribed by the functional 



TOs,[p L ]) =F[p L ,p s ] -p S / PS, Ms 
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Pl(S) being the density profile of the large (small) par- 
ticles, and F the Helmholtz free-energy functional. The 
equilibrium equation for ps in (Q) is usually a nonlinear 
equation involving both densities; the infinite asymmetry 
limit then allows to explicitly obtain p$ as an expansion 
with terms depending on pi,. This operation transforms 
(|l|) into an effective Helmholtz free-energy functional for 
the depleted large-component fluid (up to some divergent 
terms which do not affect the phase behaviour), with the 
fugacity of the solvent playing the role of an inverse tem- 
perature. 

The direct correlation function (DCF) of the effective 
fluid is defined as minus the second functional derivative 
of the excess part of ([!]) with respect to pl- ft is straight- 
forward to see that this leads to the following relationship 
with the components of the DCF of the mixture, written 
in Fourier space: 



C eS (k) = C LL (k) + 



PsC LS (k) 2 
l-p S Css(k)' 



(2) 



In terms of the structure factor this expression is equiva- 
lent to assuming Scs(k) — Sr,L{k). This is precisely the 
procedure used in |13|| to map a binary HS mixture onto 
an adhesive HS fluid in the limit e = <js/o~l — ► 0, <Jl(S) 
denoting the diameter of the large (small) spheres. 

To proceed with the functional mapping we have de- 
scribed above we need a free-energy functional able to 
deal with strong inhomogeneities (the solvent can be re- 
garded as a fluid confined in between the large particles). 
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Fundamental measure functionals are perfect candidates 
||l^ - |l7f b oth because of their ability to reduce dimension- 
ality and because they are formulated directly for 
mixtures. On the other hand, for the sake of simplic- 
ity we have chosen to apply this procedure to a binary 
mixture of parallel hard cubes (PHC). This fluid has also 
been shown to have a spinodal instability both in sim- 
ulations on a lattice p8| and in continuum space [fl6|| , 
thus its phase behaviour can be expected to be qualita- 
tively similar to that of HS. The procedure applied to 
the latter, however, has a higher level of complexity both 
because of the more complicated shape of the function- 
als jl5j and because the weaker depletion in this fluid 
requires to expand to a higher order in e. There is noth- 
ing fundamental, though, which prevents its application 
to HS, and we are currently working along this line. 

In order to determine the appropriate scaling of the 
solvent packing fraction rjs we first take the e — > limit 
in (^) forcing it to produce a finite nontrivial result. The 
right choice is 775 = O(e); thus, if we define 775 = e£, 
the resulting DCF will be (ctl = 1 in what follows) 
C ff = Cphc + C ad , where 



C a d(r) 



2(1-7?) 



{S coat (r)+yS(v)+6y 2 V(v)}; 



(3) 



CpHc( r ) is the DCF of a one-component PHC fluid 
jl6],[l7| of packing fraction 77 = t]l (actually the total 
packing fraction in the limiting fluid); z = e 2 exp(/3/j cx ) = 
£/(l — 77) is the "renormalized" excess (over ideal) fugac- 
ity of the solvent; y = 17/ (1 — 77); and 
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V(r) 


= L(x)L(y)L(z) , 


(6) 


with L(u) 


= (1 - M)©(1 - M) and e 


\v) = 1 if v > 



and = otherwise. 

Apart from corrections proportional to the overlap sur- 
face (||) and volume (||), at larger densities, the lowest 
density correction (Q) is a. delta on the contact surface — 
also proportional to the actual contact surface. This has 
its origin in the effective adhesive potential on the sur- 
faces of the cubes induced by the solvent, a nd it turns 

We 
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out to be a natural extension of that for HS 13 
will thus refer to this fluid as the parallel adhesive hard 
cube (PAHC) fluid. 

Once we know the scaling of the density we can deter- 
mine the free energy with the procedure described above. 
The fundamental measure functional of a mixture of PHC 
is given by (3F = J dr {$ id (r) + $ cx (r)}, where 



$ id (r)=^p i (r)(log[V i A(r)]-l), 



$ ox (r) = -n ln(l - n 3 ) 



ni • n 2 , n 2 , x n2, v n 2 . 
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Vi being the thermal volume of species i, Pi(r) its 
local number density, and n , n 1; n 2 , n 3 are a set 
of "fundamental" weighted densities generically de- 
fined as the convolution n a = pi * uxf 1 ' , where 
= «5ftf*f, wj x) = (6f5V5?,6!6V6?,5!5V6!), w 2 (2) = 

(W?,w?,wn> ^ (3) = nnn, with n = 

6(^/2 - \u\), and 5f = (1/2)5(^/2 - |«|). 

According to ([!]) the equilibrium scaled density of the 
solvent £(r) = e 2 ps(r) is given by 



log£(r)=/3/4 x + lne 2 - 



4 a) (r). (9) 



^ dn c 



If we now take into account that, for any function /(r), 
/ * u>^ — e a f + 0(e a+2 ), for a — 3 or any vector com- 
ponent of a = 2, 1, and /*4 0) = /+ (e 2 /8)V 2 / + (9(e 4 ), 
then £(r) = z[l -n 3 (r)] + efi(r) + e 2 &(r) + 0(e 3 ) is the 
e-expansion of ^ in terms of the n a and z (for the sake of 
brevity we omit the expressions of £1 and £2)- By over- 
line we are denoting weighted densities of only the large 
component. 

We can now insert the expansion of £ into the defini- 
tion of T [Eq. (|l|)] and expand the whole functional in 
powers of e. This leads to 

T(/xs, [pl]) = -n (e)V + /zo(e)iV + F PAH c + 0(e), (10) 
PF PAHC = [ dr {¥ id + $ cx + $ ad }. (11) 



The term —HqV + pqN diverges with e — > 0; nevertheless 
it is irrelevant for the phase behaviour because it simply 
adds Ho to the pressure and po to the chemical potential 
of the large particles. These are two divergent constants 
which just cancel out in the equilibrium equations. The 
functions $i d and <& cx are given by (0) and (§) but only 
for the large component, and <I> ad is given by 



4> 



ad 



z |Vti 3 | 2 - 4n 2 • n 2 
8 1 -n 3 ' 



(12) 



A few things are worth noticing here. First of all we 
can see that out of this procedure a new weighted den- 
sity has come up — namely Vn 3 — which was missing at 
the beginning. Secondly, the second functional deriva- 
tive of the interaction part of Fpahc yields the DCF (||), 
and this provides a consistency test. Finally, $ ad is a 
negative contribution — it arises from an attractive inter- 
action, as it can be easily realized particularizing for the 
uniform fluid (it can also be checked that this term is 
negative for any density profile). 

The equation of state of the PAHC uniform fluid, in 
terms of y = 77/(1 — 77), turns out to be 



PP = y + 3(1 - z/2)y 2 + 2y 3 



(13) 



As in the case of the adhesive hard spheres (AHS) fluid 
Q this equation has a van der Waals loop, and hence 
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yields a gas-liquid phase transition (see Fig. 1) with 
a critical point at z c = 2(1 + y/2/3) ~ 3.63, rj c — 
(1 + V6)" 1 w0.29. 
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FIG. 1. Solvent fugacity, z, vs. solute packing fraction, rj, of 
the infinitely asymmetric binary mixture of PHC, both with- 
out (a) and with (b) polydispersity (Aa/a — 4.5%). (a) The 
thick solid line separates the unstable region (U) from the 
metastable one; the thin one marks the (continuous) transi- 
tion from a metastable fluid (MF) to a metastable solid (MS); 
the dashed one is the fluid-fluid spinodal. (b) The thick solid 
line marks the fluid-solid or solid-solid coexistence; the thin 
one marks again the (continuous) fluid-solid transition be- 
low the coexistence region; the dotted line is the metastable 
fluid-fluid coexistence. 

Allowing for inhomogeneous phases, however, changes 
dramatically the phase behaviour. As usual in density 
functional theories we determine the solid free energy 
by writing pl{?) as a sum of gaussians centered at the 
nodes of a (simple cubic) lattice, leaving their width as 
a variational parameter. We also account for vacancies 
by multiplying by an occupancy probability v, which 
acts as a second variational parameter. This occupancy 
can be expressed in terms of the lattice parameter, d, 
as v = rjd 3 ; hence d can be treated as an alternative 
variational parameter — subject to the constraint d > 1. 
First thing we can see is that the fluid-solid transition 
is continuous therefore we actually obtain it directly 
from the DCF (||) by a standard bifurcation analysis — 
and that it crosses the fluid-fluid coexistence at a packing 
fraction r\ < rj c . Hence the fluid- fluid transition is pre- 
empted by freezing. 

But having d as a variational parameter reveals that 
any of these phases is metastable. The reason is that at 
any value of z, <I> a d — * — oo as d — > 1. Thus the only ther- 
modynamically stable phases this system possesses are a 
closed-packed solid coexisting with a zero-density fluid. 
This is not surprising if one compares with the phase be- 
haviour of adhesive hard spheres (AHS). Computer simu- 
lations of a square- well fluid with well-range going to zero 
|To[ reveal that the AHS fluid also has the same phase 
behaviour. This singularity was already noticed by Stell 
pot , who proved that the partition function of the AHS 



fluid diverges for a system of 12 or more particles (12 
is the coordination number of a fee lattice). It simply 
expresses the fact that at any density and adhesiveness 
(solvent fugacity in our case) the system collapses into a 
close packed solid lattice. 

Nevertheless, a more detailed exploration of the phase 
diagram shows that the fluid or the normal solid are 
sometimes local minima of the free-energy functional, 
and thus they are metastable phases, mostly separated 
from the close-packed solid by large free-energy barri- 
ers. The limit of mechanical stability of these metastable 
phases is shown in Fig. 1(a). It can be seen that the 
larger z the larger the density range of metastability. 

Thus the situation one will typically find in the PAHC 
fluid is the following. At large values of z the fluid will 
quickly collapse into the close-packed configuration. At 
small values of z, however, the fluid will be trapped in 
metastable phases for very long times; upon increasing 
r\ the fluid will be seen to undergo a continuous fluid- 
solid transition at r\ « 0.3, and it can remain solid for a 
very long period of time; if we move at still higher rj the 
system will eventually collapse, but as it may take long 
for the whole system to do it there may be an apparent 
expanded solid-collapsed solid coexistence. This picture 
is a caricature of the typical phase diagram of a fluid 
exhibiting an isostructural solid-solid transition Jl^,|ll| . 

A similar behaviour has been observed for the AHS 
fluid pl[ : the phase diagram reveals a freezing transition 
in spite that the free energy becomes concave beyond 
a certain packing fraction (depending on the adhesive- 
ness) . The latter was interpreted in j2l) as a percolation 
transition; however it is simply a sign of the collapse sin- 
gularity, which did not show up because the density of 
vacancies was fixed to (and hence the lattice parameter 
is forced to be d = r/ -1 / 3 > 1; we observe the same be- 
haviour in our system when we fix v = 1). The reported 
phase diagram is thus metastable (compare the qualita- 
tive agreement with the one shown in Fig. 1(a) for the 
PAHC fluid). 

The collapse of the adhesive potential can be removed 
by adding a small amount of polydispersity (2^j2^] . To 
see what the phase diagram looks like we have intro- 
duced a small amount of size polydispersity in the large 
cubes. In the fundamental measure formalism this simply 
amounts to replacing the weights oA by their averages 

uj 1 ^ over the chosen size distribution. We have made the 
simplest computational choice for the latter: we have re- 
placed the large cubes by parallelepipeds where each of 
their edge-lengths is chosen random and independently 
from a Gaussian distribution of mean 1 and deviation 
Act. It turns out that this choice transforms the former 
weights into "smoothed" versions of them, and that this 
removes the divergence of $ a( j no matter how small Act 
be. Besides, this choice leaves the free energy of the uni- 
form fluid unchanged. 

The results are plotted in Fig. 1(b). It can be seen that 
now we really have the typical phase diagram of a fluid 
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exhibiting an isostructural solid-solid transition — except 
that for this system the fluid-solid transition is continu- 
ous. For large values of z we have a very diluted fluid 
and a very dense solid separated by a large coexistence 
region. As we decrease z the coexisting fluid eventually 
becomes denser until it reaches the density at which freez- 
ing occurs. From then on and down to a critical value, 
z s , coexistence is between an expanded and a dense solid 
with the same crystal structure. 

The results we have obtained can be reinterpreted in 
terms of the original very asymmetric binary mixture. 
The fluid phase would correspond to the stable mixture, 
and the solid phases to the solid of large particles perme- 
ated by a fluid of small particles. According to the PAHC 
phase diagram — either mono or polydisperse — no fluid- 
fluid phase separation will ever occur for very asymmetric 
binary mixtures of PHC; instead, at large solvent fugac- 
ity (large solvent density) there will be fluid-solid phase 
separation, and at smaller solvent fugacity the phase sep- 
aration will take place after the large particles have crys- 
tallized. 

Notice that the existence of this fluid-solid phase sep- 
aration largely increases the region of the phase diagram 
where the mixture is unstable. This explains why recent 
simulations on very asymmetric binary mixtures of PHC 
H report the mixture to be unstable at packing fractions 
of the large component much smaller than those pre- 
dicted by fluid- fluid phase separation jl(| (notice that, 
according to the authors, the nature of the coexisting 
phases cannot be discerned in their simulation). 

The same kind of phase behaviour has been very re- 
cently reported to occur in simulations of binary mixtures 
of HS at finite (but small) diameter ratios using the deple- 
tion potential (they also show a solid-solid transition 
in a certain range of solvent fugacities). We can see then 
that this phase behaviour can be inferred from that of the 
AHS fluid [|T| . This result is of great importance because 
it establishes a connection between the phase behaviour 
of mixtures and that of the solvent fluid interacting via 
an adhesive- like potential |lO|,[Ll]], thus opening a new 
route to study demixing in hard-particle fluids. This re- 
markable enhancement of the slow divergence of the pair 
correlation function of two large bodies at contact, in- 
duced by depletion jl2| , is probably due to the fact that 
the confinement of the solvent in between the solute par- 
ticles induces strong inhomogeneities in the former [see 
Eq. (^)] an hence largely increase depletion. 
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